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WHAT IS CLAIMED IS : 

1. A method for the fourth-order, blind identification of at least two sources 
in a system comprising a number of sources P and a number N of reception 

5 sensors receiving the observations, said sources having different tri-spectra, 
wherein the method comprises at least the following steps: 
a) a step for the fourth-order whitening of the observations received on the 
reception sensors in order to orthonormalize the direction vectors of the 
sources in the matrices of quadricovariance of the observations used, 
10 b) a step for the joint diagonalizing of several whitened matrices of 
quadricovariance (step a) in order to identify the spatial signatures of the 
sources. 

2. A method according to claim 1, wherein the observations used 
15 correspond to the time-domain averaged matrices of quadricovariance 

defined by: 

p 

£>v(T|,T2,T3) = ^(ti,T2,T 3 ) {a p ® a/) (a p ® a/) H (4a) 
p=i 

= A Q C 5 (t,,t 2 ,t 3 ) A q h (4b) 
where A Q is a matrix with a dimension (N*x F) defined by A Q = [(ai®ar), 
20 (a p ®a p % Cs(ti,t2,T3) is a diagonal matrix with a dimension (Px F) defined by 
Cs(xi,T2,t 3 ) = diag[ci(Ti,T 2 ,T 3 ),..., c p (ti,t 2 ,t 3 )] and c p (xi,T2,T3) is defined by: 

Cp(Ti,T2,T 3 ) = <Cum(jp(/), 5p(M,)*, Sp(t-T 2 )\ Sp(t-T 3 ))> (5) 

3. A method according to claim 2, comprises at least the following steps: 

25 Step 1: the estimation, through Q;\ of the matrix Q Xi from the L 
observations x(lT e ) using a non-skewed and asymptotically consistent 
estimator. 

Step 2: the eigen-element decomposition of Q; A X , the estimation of the 
number of sources P and the limiting of the eigen-element decomposition to 
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the P main components: Q A X « £;\ A; A * £; A * H , where A; A X is the diagonal 

matrix containing the P eigenvalues with the highest modulus and E A X is 

the matrix containing the associated eigenvectors. 

Step 3: the building of the whitening matrix: T A = (A; A x)~ I/2 E A ^. 

Step 4: the selection of K triplets of delays (xi k y T 2 k ,T3 k ) where 

|ti*| + |t 2 *| + |t3M*0. 

Step 5: the estimation, through Q; A x (xi k , 12*13*), of the K matrices 
<?x(t 1*^2*13*). 

Step 6: the computation of the matrices T A Q; A x(Ti k ,x 2 k } X3 k ) T A " and the 
estimation, by U' A so i, of the unitary matrix U so \ by the joint diagonalizing of 
the K matrices T A Q; a *(ti* T 2 fc ,T 3 fc ) T A " 

Step 7: the computation of T A V; A S0 <=[b; A i.. .£>;%] and the building of the 
matrices B; A , sized (Nx TV). 

Step 8: the estimation, through_a; A F , of the signatures a, (1< q < F) of the P 
sources in applying a decomposition into elements on each matrix B\ A i 

4. A method according to claim 1 to 3, comprising at least one step for the 
evaluation of the quality of the identification of the associated direction 
vector in using a criterion such as: 

D(A, A) = (oti, ct2, ,<xp) (16) 

where 

a p = ™™l d ( a p> 4)1 (17) 

and where d(u,v) is the pseudo-distance between the vectors u and v, such 
that: 

I H I 2 

d(M ' = 1 - riroo (is) 



28 



5. A method according to claim 1, comprising at least one step of 

fourth -order cyclical after the step a) of fourth-order whitening. 

6. A method according to claim 5, wherein the identification step is 

performed in using fourth-order statistics. 

7. A method according to claim 1 wherein the number of sources P is 
greater than or equal to the number of sensors. 

8. A method according to claim 1, comprising at least one step of 
goniometry using the identified signature of the sources. 

9. A method according to claim 1, comprising at least one step of spatial 
filtering after the identified signature of the sources. 

10. A use of the method according to claim 1 in a communications 
network. 



